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ABSTRACT 


Fluctuations of the East-West component of the Earth's 
geomagnetic field were measured at a remote land site. The 
resulting data were transmitted by a solar powered telemetry 
system to the Naval Postgraduate School at Monterey, California. 
BHoeche pOWer spectra for the frequency range of .02 - 20 Hz 
calculated. The measurements, which covered a 4-month 
meecrvyal (July 20 - October 10, 1981), consistently show a 
menimum Of Activity in the interval 3 - 7 Hz. At frequencies 
below the minimum, in the range of .02 - 3 Hz, the typical 
monotonic decrease in background activity with frequency was 
observed. At 1 Hz an average power spectral density of 1X107? 
nT’/Hz was observed dumm the day and 3.1X10>° n?*/Hz at night. 
In contrast, at frequencies above the minimum, in the range 
fee Hz, the activity is dominated by the first Schumann 
resonance. An evaluation of the East-West component spectra 
and concurrent underwater horizontal component measurements 
showed a 90% correlation with the underwater spectra. The 
underwater field strengths were normally 3-5 dB less than the 
Strengths measured on land. Splitting of the first Schumann 
meoomanee peak Into a doublet structure was observed in 10% 


Bf the Land data. 
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I. INTRODUCTION 


This thesis research is part of an ongoing effort at the 
Naval Postgraduate School to obtain improved long-term data 
and interpretations of the electromagnetic noise on the ocean 
floor. The project objective is to study and interpret 
Peotlais im the frequency range from .01 Hz to 100 Hz. The 
overall project emphasizes the importance of obtaining measure- 
Memes Of geomagnetic noise fluctuations on the sea floor and 
on land over a period of several years. 
The particular objectives of this thesis are: 
1. To install a remote land-based monitoring site and 
a telemetry system for the transmission of geomagnetic 
data to the Naval Postgraduate School. 
ee OeGOMlCCENthe CaSt=yoummeomponent of the local 
Morte rela T1UCceiaeion data for comparison with 


sea floor data and to test the system's sensitivity. 





II. BACKGROUND 
A. MAIN MAGNETIC FIELD 


Lom@memmapmetlGsolrveys at ground level and, more recently, 
by satellite indicate that the principal source of the main 
Beomagnetic field is beneath the earth's crust. Approximately 
90% of the main field exhibits characteristics roughly 
equivalent to the field of a terrestrial-centered short bar 


magnet or dipole inclined at 11.5% to the axis of rotation. 





Figure 1. Dipole Appearance on Geomagnetic Field 
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Peeves OF THE FIELD 


MicmvangnOlsmeomponemes Of the magnetic field are shown in 


Figure 2. They are defined as: 
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North-South Component 
East-West Component 
Vertical Component 
Declination Angle 
Hom@izontal Component 
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Inclination or Dip Angle 


ZENITH (-) 


sien “. NORTH (+) 





Fisure sc.) Elenie@mts of the Field 
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Pee mscnOreticmramiliar Dar magnet, a magnet perfectly 
meee tomturm in any direction would show the direction of the 
melas oince 1t is almost impossible to suspend a single 
magnet in this manner, two are used: one with a vertical 
pease (cOmpass) and one with a horizontal axis (dip needle). 
The compass shows the direction of the horizontal component H; 
and when properly oriented, the dip needle shows the inclina- 


Pee Or dip. 


eet iME VARIATIONS OF THE FIELD 
The time variations of the geomagnetic field can be 
meeadly Categorized into “quiet variation fields", "disturbed 
fometation fields’, and “micropulsations"'. 
MOLEC cmVar tactom Field 
Smoothed geomagnetic records show a consistent trend 
that clearly indicates a pattern of daily variation with 
meepect to solar local time. The average variation patterns 
derived from chosen quiet day records define a variation 
meeld Called the SOLAR QUIET DAILY VARIATION which is denoted 
Sage there is also a lunar daily variation, denoted L, but 
1t is much weaker (L<0.1 Sq.) and more variable. 
ieGecthtirds OL these daily (diurnal) fluctuations, 
me., ehose having a 24 hour period, are due to current 
sources external to the earth. The remaining one third is 
due to internal currents induced in the earth's surface 


layers by variable external fields. The external currents 


LZ 


responsible are found to flow at an approximate altitude of 
100 kilometers, and are produced by convective movement of 
charged particles moving across the earth's magnetic field 
lines. The convective motion of the particles is caused 
mainly by solar heating of the upper atmosphere. 

Therefore, the sun produces a stationary current 
system in the upper atmosphere, and the earth rotates under 
it once a day. This mechanism, known as the Atmospheric 
Dynamo, starts in Monterey at about sunrise and reaches its 
peak about noon. Sq. varies characteristically with latitude. 
At the equator the maximum horizontal intensity is 100 nT’, 
where at higher latitudes it is generally -25 to -50 nT. 

Sq. is also dependent on the season and the phase of the 
solar cycle. In the summer and during a sunspot maximum, 
Sq. 1S increased; in the winter, it is decreased. 

2. Disturbed Variation Fields 

Any fluctuations other than the quiet day solar and 
lunar variations are called magnetic disturbances. All 
disturbances, other than those directly attributable to a 
solar flare or fluctuations produced by upper atmospheric 
irregular motions, are caused originally by disturbed solar 
plasma and may therefore be generally classified as belonging 
to one family. 

[Irregular magnetic fluctuations of one-to-two hours 
duration are related to changes in the solar wind pressure on 


1 - -5 
imal: 10 Gauss 13 





magnetosphere and associated changes in the southward component 
of the interplanetary field. These fluctuations are sudden 
impulses and magnetic bays. ‘Sudden Impulses' are sudden 
increases of several nanoteslas followed by a gradual increase 
or decrease in the field. This is succeeded by a return, 

With some small oscillations, to the normal field. These 
latter fluctuations usually last about one to two hours. 

The term magnetic storm is reserved to apply to a 
relatively severe, long lasting disturbance with recognizable 
features. Magnetic storms are caused by plasma bursts from 
the sun, since a solar flare emits both X rays and plasma. 
time X rays precede the plasma and enter the earth's atmosphere, 
causing increased ionization in the sunlit atmosphere. This 
ionization occurs in the lower ionosphere and enhances the 
currents in the atmospheric dynamo thereby producing what is 
known as the Solar Flare Effect, a precursor to magnetic storms. 

The major magnetic storm disturbances are caused by 
dynamic pressure changes of the solar wind. When a plasma 
blast from the sun arrives, it suddenly increases the solar 
Wind pressure, compressing the magnetic field (Sudden commence- 
ment). It maintains the compression for a time, called the 
initial phase. When the pressure decreases with the passing 
of the plasma blast, the circulating currents which have been 
induced in the magnetosphere tend to make the field intensity 
"overshoot" before the currents gradually dissipate. This 


initiates the recovery phase back to the normal quiet field 
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condition. The entire magnetic storm can last from one to 
three days with a recovery time that may be even longer. 
wemocomagnerae Micropulsation 

Field changes which occur at periods of 0.2 seconds 
to 10 minutes are called 4icropulsations. These characteristi- 
cally have amplitudes from tens of nanoteslas to a fraction 
mamone Manotesla. Micropulsations are generally divided into 
two types. The terminology recommended by the International 
Association of Geomagnetism and Aeronomy is pce for pulsation 
Seimemnuous and pi for pulsation irregular. The continuous 
(pc) micropulsations have amplitude variations which are 
quasi-sinusoidal. The irregular (pi) micropulsations exhibit 
/meeetitaritzes in both frequency and amplitude. The period 


Pemses and average amplitude for these nulsations are: 


pet Boome OP 
pe2 5-10 sec Ook Lon T 
peo 10-45 sec GAte aa nt 
pc4 45-150 sec ee 
pcs 150-600 sec Poe Ont 

Sail 1-40 sec Pie a 0s bm, 
joes 40-150 sec dees) et 


Micropulsations below 3 Hz are produced mainly from wave- 
particle interactions in the magnetosphere. Pulsations of 
meequencies from 5 to 3000 Hz make up the ELF region of 
meeere 5. This ELF region has three principle contributors: 
mee orerics’, ELF Emissions, and Earth-Ionosphere Cavity 


Resonances. 
15 
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Pepmomerites "ame Clectromagnetic signals from atmospheric 
electric discharges that propagate in the earth's waveguide 
between the ground and the lower boundary of the ionospheric 
E-region. The "Sferic' wave travels very long distances and 
its waveform consists of a main high frequency oscillatory 
head (VLF), followed by a lower frequency (ELF) tail-like 
oscillation, sometimes referred to as a "slow tail". "Sferics" 
commonly last about 20 milliseconds and have frequency 
components from 30 to several Hertz. 

ELF emissions are an excitation of whistler waves 
(300 to 30,000 Hz), made by charged particles streaming along 
the earth's field lines. The whistlers in the VLF region 
sometimes produce lower frequency components in the ELF 
region. 

The Cavity Resonance signals are resonantly excited 
Pee tehtnineg transients in the concentric spherical cavity 
between the earth's surface and the lower region of the 
ionosphere. The power spectra of the signals show maxima 


ete ao. 14.1. 20.3, 26.4 and 32.5 Hertz. 


tee PREVIOUS WORK 

PRemiOUsewerlh by Santirocco and Parker (1965) and Park 
(1964) indicates that the geomagnetic field decreases with a 
Slope of approximately -6 dB/octave between 10°*and 1 Hz. 
mom i to 40 Hz, the field exhibits a levelling off and is 


dominated by the Schumann Resonances (Schumann and Konig, 1954). 
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Recent measurements in the frequency range of 0.1-14 Hz 
mememeomauerca by A, €. Fraser-Smith and J. L. Buxton (1975). 
Their measurements were taken over a two month interval at 
Stanford, California which is 100 km from the site used in 
this work. They found that the general decline in the slope 
(-6 dB/octave) continued to approximately 5 Hz, where the 
decline was arrested by the superimposed Schumann resonance 
activity. Their results indicated the first Schumann resonance 
aeeapproximately 8 Hz. 

More recent work was accomplished in the low frequency 
ranges at the Naval Postgraduate School in Monterey, California. 
Barry (1978), Clayton (1979) and McDevitt and Homan (1980) 
investigated geomagnetic fluctuation in the frequency ranges 
memuelz, .4-40 Hz and .04-25 Hz, respectively. Their results 


@eenwin general agreement with earlier measurements. 
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Cewes ocm  DEoCRIPTION 


Eee SEMOTE MONITORING SITE 

The remote geomagnetic monitoring station was established 
primarily in support of the Naval Postgraduate School's 
research work on geomagnetic fluctuations on the ocean floor. 
To measure magnetic fluctuations on the sea floor, a source 
of reliable, nearby land reference data is needed. Precision 
time and frequency correlation of sea floor data with land 
itweaepermits differentiation between the local sea-generated 
fields and those due to external sources. The remote monitor- 
ing site is located in the mountain 40 km. southeast of Monterey, 
California, in the Carmel Valley near the Tassajara Road 
entrance to the Los Padres National Forest (Figure 4). The 
hilltop altitude of 3844 feet with the additional 20 feet of 
antenna provides the needed elevation for VHF (138.7 and 140.2 
MHz) propagation to the Naval Postgraduate School. This 
Mecation is a Magnetically quiet site, relatively close to 
the ocean (25 km.). The lack of low frequency industrial 
and power transmission interference coupled with an abundance 
of solar energy for power were also considered in the site 
selection. The basic telemetry design for the monitoring 
station was performed by Beliveau (1980). 

At the base of the antenna, a storage cabinet was constructed 


to house the radio and electrical equipment. A solar array 
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was mounted on top of the cabinet and oriented for maximum 
power output. A single coil sensor was buried in a wooded 
area on the far side of the hillcrest, approximately 90 


Mmemwers from the transmitter location. 


Eee oYSlEM COMPONENTS 

RiemimlajyOb Parts Of the system are shown in Figure 5, 6, 
and 7 and further details are given in Appendix A. 

The sensor, a coil antenna, was oriented to measure 
fluctuations in the east-west component of the local magnetic 
field. As the field varies, small voltages are induced in 
the coil. These voltages are amplified and transmitted by 
cable to the storage cabinet where they are amplified a second 
time to the desired signal level. The resultant voltage 
variations are applied to a voltage controlled oscillator. 

As the magnetic field changes, it causes the VCO to vary 
about its free running frequency ( 1500 Hz). This time 
varying frequency is then transmitted by radio to an antenna 
atop Spanagel Hall at the Naval Postgraduate School. 

The incoming signal is detected with a phase-lock loop 
that tracks the tone and removes any static bursts. The 
phase-lock loop frequency is recorded for later analysis. 

During analysis the tape recorder output is demodulated 
with a frequency to voltage converter (FVC) producing a 
varying voltage directly proportional to the original output 


from the sensor coil. The VCO/FVC combination introduces a 
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a positive dc voltage offset that is removed with a differen- 
tial amplifier. The resultant a.c. signal is processed by a 
spectrum analyzer that produces a 400 data point digital 
display. Data points thus obtained are entered into a Hewlett- 
Packard HP-85 computer. The HP-85 corrects each data point 

for the overall system gain function as well as for the 
e@ectrum analyzer bandwidth. (Section III, C and D). The 
mamal processed data is plotted with a Hewlett-Packard /225A 
Meoturer. A more detailed discussion of the pieces of equip- 


ment and their characteristics is contained in Appendix A. 


Mees Ol EM TRANSFER FUNCTION 

The entire system can be represented by a block diagram 
that converts a fluctuating magnetic field into a voltage 
that can be displayed and whose power spectral density can 


be measured with a spectrum analyzer. 


B(f) SSieM Kaeo) 


(nanoteslas) Ga (volts) 





Since the transfer function, H(f), is frequency dependent, 
it was calculated by dividing the system into components 
as shown below. The sensitivity of each component was found 
over the frequency range of interest and then added point by 
point to produce the system transfer function. 

The radio, phase-lock loop, tape recorder and differential 
amplifier are all of unity gain and do not affect the transfer 


function calculation. The data presented in this research is 
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nonereneceue tO LF tl in order tO be consistent with earlier 
work. Details of the procedure used are described in 


Appendix B, Sections I and II. 
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Given the system transfer function, the corrected power 


spectral density can be found from 
P = N - 10 log [H(f)]? 


where 


N = The Spectrum Analyzer reading in db, referenced to 
i Vole 


HW = The Hanning window bandwidth correction for the 
particular Spectrum Analyzer used in db 


Fee ) The transfer function defined above 


P = The corrected power spectral density where 0 db is 
referenced to 1 nT*/Hz 

ne SYSTEM NOISE 

The major source of the total system noise is that intro- 
duced by the preamplifier located at the sensor. The best 
amplifier performance is achieved when the amplifier causes 
the least decrease in the overall signal-to-noise ratio. The 
preamplifier in this system, Model 13-10A, was selected because 
it possesses a better signal-to-noise ratio and improved low 
noise characteristics when compared with Model 15-10 which 


was used in the earlier work of McDevitt and Homan, (1980). 
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To compare the two preamplifiers, a 120 2 resistor was 
connected across the input while the rms output voltage was 
measured in the spectrum analyzer. Figure 9 shows that 
approximately 10 dB of improvement was obtained at .1 Hz by 
using the Model 13-10A. 

The total system noise rms voltage level was measured 
several times. A120 2 resistor (to replace the 124 2 sensor 
coil) was connected to the preamplifier input and normal 
analysis of the resulting transmitted signal was performed. 
The resulting system noise obtained is shown in Figure 10. 

A comparison with the noise measured by McDevitt and Homan: 
at Chew's Ridge shows that the telemetry system noise was 

15 dB + 5 db lower over the frequency range of 1-20 Hz. As 
discussed above the main reason for this difference is the 
use of the improved preamplifier. The remaining improvement, 
reduction in system noise, is due mainly to the substitution 
of an instrumentation tape recorder for the portable stereo 


cassette type used by McDevitt and Homan. 
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ieee eR IMENTAL RESULTS 


A. INTRODUCTION 

During the months of July, August, September, and October, 
mie, a total of 68 hours of data was collected. Four specific 
mememintervals, 0730-0900, 1100-1230, 1800-1930 and 2300-0030 
were selected each day. Approximately the first 45 hours of 
data were subject to static bursts in the telemetry transmission 
that prevented analysis by long (N greater than 8) time 
averages of the data. During this initial installation phase, 
numerous equipment adjustments and modifications were made. 

The data presented represent measurements obtained from the 

final system (4-5 October, 1980) and data taken concurrently 
with undersea measurements made by (Ames and Vehlsage, 1981) 
cums and 24-25 August, 1981. 

Comparisons were made of the power spectral density for 
each of the above time frames for a number of days. Comparisons 
were also made for spectra at different times of the day to 
determine daily variations. The data was correlated with 
previous data to show compatability for each time period 
using quiet days selected from the data. The days of varying 
degrees of magnetic activity were investigated separately. 

The level of magnetic activity was determined by using 


the daily 24-hour A-index and the 53-hour K-indices measured 
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at the Fredericksburg, Virginia (mid-latitude) magnetic 
[Pomerorimneestatton. The Space Environment Services Center 
located in Boulder, Colorado, publishes a Preliminary Report 
and Forecast of Solar Geophysical Data, including geomagnetic 
activity. A more detailed synopsis of geophysical data is 
contained in the environmental data services montlhy publication 
"Solar-Geophysical Data". The A-indices range from 0 (very 
quiet) to 400 (extremely disturbed) with an A-index of 30 
or greater, indicating geomagnetic storm conditions. The 
K-indices range from 0 (very quiet) to 9 (extremely disturbed). 
The position of the Chew's Ridge sensor was fixed, to 
measure the East-West component of the local magnetic field, 


geesnown in Figure il. 


B. REVIEW OF DATA 

Figures 12-27 represent eight separate data tapes. Each 
data tape was analyzed over two frequency scales, .01-5 Hz 
apae.t-20 Hz. Due to differing amounts of time required to 
produce each average, a larger portion of the tape was analyzed 
for the .01-5 Hz average. In an effort to minimize differences 
between the two scales, the .1-20 Hz average was made from 
the approximate center of the tape used for the .01-5 Hz 
average. Figure 28 is an average of all daytime and nighttime 
data recorded during the project. 

The data presented in Figures 12-19 represent four 


measurements taken concurrently with undersea measurements 
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Figure 11. Chew's Ridge Sensor Orientation 
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made by (Ames and Vehlsage) on 21 and 24-25 August from the 
Research Vessel Acania. The nature of the underwater data 
collection system made it necessary to deviate from the 
normal collection schedule. The underwater system, when 
placed on the sea floor, used a timer to start its recorder 
after the ship had moved far enough away to no longer affect 
the readings. Once a thirty-minute tape had been recorded, 
the system had to be recovered, reset and once again lowered 
to the ocean floor. The normal two hour turnaround time for 
@mmsOperation was Subject to numerous possible delays. In 
order to make the concurrent measurements, the Research Vessel 
Acania, radioed each start time to the Chew's Ridge recording 
Site in Spanagel Hall at the Naval Postgraduate School. The 
land recordings were started within + 30 seconds of the sea 
recordings, in all cases. 

The August measurements were made prior to the addition 
of a 14 dB amplifier that was included in the final system; 
therefore, the transfer function used in processing this 
data was 14 db lower than the one given in Section III. The 
power spectral density for each E-W Chew's Ridge measurement 
was compared to the results of the horizontal component 
measurement obtained from the undersea measurement. 

Three of the four land-sea measurements were in excellent 
merecmeme. Figures 13, 17 and 19 (the midday, night and 
Sumrise spectra for the .1-20 Hz scale) are consistently 3-7 


db higher than those measured underwater. In addition 
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individual peaks observed exhibit a 90% correlation. A differ- 
ence in relative magnitude between the peak in the two sets 

of data was noted. This is probably due to the fact that 

the land data was fixed to measure the East-West component, 
while the sea data system measured a horizontal component 

of unknown orientation. 

The fourth land-sea comparison, Figures 14 and 15, is 
interesting in two aspects. First, the power spectral density 
on the .1-20 Hz scale (Figure 15) for the local sunset, 1935 
hours on 24 August, is the lowest of the day. On the .01-5 Hz 
scale however, the measurement crosses over the other measure- 
fetes rOr the day for frequencies less than .1 Hz and 
reaches the highest level measured to date on the system. 
Secondly, the fact that this measurement occured on a day 
with an A index of 15 and K's of 3's made this spectrum suspect. 
The concurrent underwater measurements also showed an even 
ereater surge of power spectral density increase across the 
entire frequency scale that was 30 dB greater than that 
recorded on the land system. The correlation between the 
individual peaks of the two sets of spectra was very poor 
and makes this result questionable. Another indication that 
this time period was nontypical is that the local night 
measurement, Figures 16 and 17, were among the highest 
recorded for the entire day, a result just opposite of what was 


normally observed. 
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niemorpeeendmitectoures Z2Z0-2/ were obtained from data 
ar ewedunane normal recording Operations on 4-5 October, 1981. 
The geomagnetic activity level was active to near storm 
levels, A = 24, and K's were mostly 4's and 5's, during 
fais period. This set of spectra is characterized by more 
prominent peaks (up to 10 dB) than was observed in the 
maemst data, A sharp peak is seen at approximately 1.8 Hz 
mimeoocn the .1-20 Hz and .01-5 Hz graphs for 0830 hours on 
5 October (Figures 26 and 27). This appears to be a large 
scale Pcl event associated with the magnetic storm activity. 
It 1s also near the typical 2 Hz peak seen in the average 
nighttime data (Figure 28). A second interesting feature in 
mitemset Of Spectra is the splitting, or apparent splitting, 
of the 8 Hz Schumann resonance peak into a doublet with one 
EU@@erarery peak near $8 Hz and the other near 10 Hz in Figure 
25 (4 October, 2330 Hours). While this phenomenon was not 
as obvious as that noted by Fraser-Smith (1975) in his North- 
south component research, it does seem to occur approximately 
10¢ of the time and to remain for several hours during each 
Securrence. 

Figure 28 shows the average East-West component during 
local day and night measurements for the period August-October, 
1981. (Local California time listed on the graphs was 7 hours 
behind UT during this time frame.) On the average there was 
a noticeable decline in the geomagnetic activity level at 


night. The decline is most often seen at the low frequencies: 


44 





(ZH S-TO°) SUOTJeNINTY platy Iryousep ysoM-qaseg 


ZH AINSNDAYS 
l T° 


re legis pla eee 


LNANOdWOO LSAM-LSVda FADGTIY S,MadHo 


SYH BETTI LI0 V 


"OZ dandy 


Td” 
BB- 


Be 


ALISNGG TWaLoddsS wdMhod 


(dP) 


ZH/2Lu 1 
>JOYN AP O 


45 





ibe. ee ent 
a re i"er I. 
_* 
ae 


| - 


(ZH OZ-1°) SUOTIeNINTY pLloty DItyousepW asom-aseq “TZ sandy 


ZH AINSNDAYS 
BT TI 


LNANOdWOO LSAM-LSWa aOdTY S,MAHO 


SYH Sell LOO 7 


ALISNGG TWwuLodds wdhOd 


(dp) 


ZH/zLU T 


-JOWY AP 0 


46 


(ZH S-10°) SUOLJENIINTY PLOTy IT IUdeW YSOM-ISeq *7Z 9INdTY 


ZH AINSNVSAYS 
I Ke 


LNANOdWOO LSAM-LSWaA ADSDGIY S,MaHO 


SYH SS8T LI0 


TA” 
BB- 


B 


BT 


Be 


ALISNGUG TWaLoddS udiiod 


(dp) 


ZH/7LU T 
:JOu AP O 


47 





(ZH OZ-1°) SuotTzZeNIINTY PLety ITyouBeyy saM-asey 


ZH AINANDAYS 
at T 


INANOdWOO LSHM-LSVd aDdTY S,MaAHO 


SYH SS8T 190 V 


"¢Z ainstyy 


AtPslaegsTVanosgds ddiOd 


(ap) 


ZH/z7LuU T 
:JOuU AP 0 


48 





(ZH S-TO°) SUOTJENIINEY PLSty DiTysusey ysom-3seq “pz asansry 


ZH 


AINANDAYS 


T 


LNANOdWOO LSAM-LSVWd dAOGIY S,MdHo 


SYH BEEC LUO PV 


Td” 
BE- 


Bt 


Be 


traboadS sxaliod 


ALISNdG 


(ap) 


ZH77Lu T 
-JOuN AP O 


49 





-_—— — ia 


(ZH OZ-1T°) SuoTJeNIINTY PLSTtY It yousep ysom-yseq  °SzZ sand, 


Bt 


ZH AINSNOSYS 


LNANOdWOO LSAM-LSVd FAOCIY S.MdaHO 


SYH BEEc LI0 


ie- 


At- 


BT 


? Ce ra 1 
eae det 


™ 
‘ 


ee ae hc 


eRe s 
i 
Bara 


, +r 
Naess 


A 


(+7P) 


ZH/7LuU 


-JeON AP 


50 


T 
0 





(ZH S-10°) SUOTJeNIINTY PLSTyY Dryoudew Asom-asey “gz sind TY 


ZH AINANDAYS 
T i 


a RN eee Wer ees en Weer cag er 0 aa 


LNANOdWOO LSAM-LSVWG dASGIY S.MdHO 


SYH BEBB LIO GS 


BT 


Be 


AL ESteG TwatlodgdsS wibiod 


(dP) 


ZH/zLuU T 


-Jeu AP O 


Dik 





(ZH OZ-1°) SUOTIENIIN{[Y PLeTY DrIRWeUseW 4saN-iseqg 


ZH AINANDAYS 
BT T 


LNANOdWOO LSAM-LSV4A AOGIY S.MaHO 


SYH SEBB 190 S 


"22 9INS EY 


ALISNGHG TWuLodds wdhOd 


(ap) 


ZH/7LU T 
Jou AP 0 


Sie 





at 0.1 Hz, the average power level drops by roughly 10 dB. 
An approximate 3 hour lag was also noted after both sunset 
and sunrise until magnetic activity moved to its normal 
nighttime and daytime levels. The average East-West spectrum 
clearly shows the minimum ULF geomagnetic activity in the 
interval 3-7 Hz as reported by Campbell (1966, 1969) and also 
by Fraser-Smith and Buxton (1975) in the North-South component. 
No other minimums were observed in the frequency interval 
0.1-20 Hz. In particular, there is no minimum in the frequency 
interval 0.2-0.5 Hz which separates Pcl geomagnetic pulsations 
meomerced (0.1-0.2 Hz) and Pce3 (0.022-0.1) pulsations. A 
minimum of activity could be anticipated near 0.2 Hz because 
of the very different properties of the Pcl and Pc3 classes 
of pulsations (Jacobs, 1970). Instead, we typically observed 
a monotonic decrease of activity with increasing frequency 
throughout the frequency range of .1-3 Hz; however, this 
could be due to the system's lack of sensitivity at the low 
frequency end of the spectrum. 

ime wtirst Schumann resonance, at approximately 8 Hz, 
produced an increase in the level of geomagnetic activity 
meerrequencies in the range of 6-10 Hz and is clearly seen 
in the average spectra. Small spectral peaks appear at 
approximately 2 and 5 Hz during the nighttime hours. These 
could mark the occurrence of small-amplitude Pcl pulsation 


events. 
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C. OBSERVATIONS AND RECOMMENDATIONS 

iteworinicmonyecs Of Lhe remote site is to provide a 
source of reliable data that can be used in support of the 
Naval Postgraduate School research work on geomagnetic fluc- 
tuations on the ocean floor. This required that the data 
produced by the system during the initial calibration period 
be checked for consistency with other sources as a measure 
Seerts reliability. 

A review of previous data taken at Chew's Ridge by Naval 
Postgraduate School students found one measurement of the 
East-West component by McDevitt and Homan in May, 1980, over 
the 1-25 Hz frequency range. Measurements made by Fraser- 
emmta and Buxton (1974) in Stamford, California, approximately 
80 miles north of the Chew's Ridge site were also reviewed. 
The Stanford data was for the North-South component, but 
was felt to be valid, basedupon work by Wertz and Campbell 
(1976). They found that the North-South and East-West field 
components were about equivalent in strength at a mid-latitude 
location ae 2-95 - A third source of mid-latitude East-West 
horizontal magnetic component measurement was Polk and 
Fitcher (1962). Each of the three sources gave an average 
ewer at wi, | and 10 Hz that was within + 15 dB from the 
average value of the East-West component obtained from the 
telemetry system. 

A difference of 3-8 dB was also noted between the value 


obtained at the same frequency when averaged over different 
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Eeequemecy scales. This effect is due to the different amounts 
of time required to obtain the same number of averages when 
using the spectrum analyzer - the smaller the frequency range, 
the greater the time required. An example of this effect 1s 
seen in Figures 16 and 17 (25 August, 1981, 0045 hours) where 
the .1 to 20 scale has a large spike (greater than 15 dB) at 
approximately 3.1 Hz. The .01 to 5 Hz scale spectra also 
Contain this spike, but it is at least than 10 dB and other 
spikes of nearly equal magnitude have appeared. By watching 
the instantaneous display on the spectrum analyzer, it was 
determined that the 3.1 Hz spike was very strong over a four 
minute period but became much weaker for the next twenty 
minutes. This indicates an inherent problem of attempting 
to use spectral analysis methods to evaluate non-stationary 
processes, such as the earth's changing magnetic field. 

While the general behavior of the observed spectra agrees 
With that illustrated in Figure 3, the amount of data collected 
during this experiment was not adequate to draw conclusions 
that could be used as a basis to evaluate underwater data. 

A longer period of continous observations, extending through- 
out the year and for different levels of magnetic activity, 
is required. A prime object of this effort should be to 
determine the optimum frequency bandwidth (averaging time) 
that can be used during different times of the day and at 


varying magnetic activity levels. 
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V. EQUIPMENT/SYSTEM IMPROVEMENTS AND RECOMMENDATIONS 


The system utilized for the project was subject to 
numerous breakdowns and was unable to meet the prime objective 
of providing continuous geomagnetic data. 

The main cause of system down time was the overheating 
of the radio transmitter. Once the maximum operating temper- 
ature was exceeded, the radio shut down automatically. 
Numerous attempts were made to reduce the radio temperature 
by shading, ventilation and attaching thermal radiators. 

When none of these attempts proved a deterrent to the over- 
heating, it was concluded that the basic equipment designed 

by Motorola was the source of the problem. The MX 320 series 
radio does not appear to have the heat sink capacity necessary 
PomallLow rapid radiation of the heat produced during continuous 
operation; therefore, overheating developes. 

Ghemneat problem 1s compounded by the fact that the 8 volt 
dc regulated Motorola solar system is designed to vary from 
8 to 12 volts dc, according to temperature. These problems 
have been reported to Motorola and an engineering solution 
has been requested. 

Until an engineering solution is received, the solar 
system could be disconnected and the radio powered by the 
storage batteries alone. This would provide a constant 


voltage of 8.0 volts or less for a pericd of two weeks at 
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which time the batteries would be 50% discharged. The 
batteries could then be recharged by transporting them back 

to the Naval Postgraduate School or by disconnecting the radio 
and connecting the solar array to the batteries for charging. 
This method would allow half time operation (assuming that 

the radio can operate at a constant 8 volts dc without 
overheating). 

For approximately $475 an additional 4 batteries could 
be purchased that would allow almost 100% operation by 
alternating the charge and discharge of the two systems from 
100% to 50%. This method would require bi-monthly maintenance 
Visits to the remote site, a significant increase in main- 
tenance from the original system design. 

Another solution would be to design a voltage regulating 
circuit between the batteries and radio that would allow only 
7.5-8.0 volts into the radio. The small difference between 
the desired voltage and the source voltage makes this difficult 
to achieve. 

Once the solar system voltage fluctuation problem is 
solved, the system could be further improved by redesigning 
aiem CO lana. 14 dB amplirier to operate on + 4 volts dc. This 
change would allow all of the electronic devices in the 
storage cabinet to be powered from the storage batteries and 
two + 6 volt Gelyte batteries could be eliminated from the 
System. It should be noted that this change would prevent 


the VCO and amplifier from being interchanged with those 
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used in the underwater experiments. At present the theoretical 
limit to the operating life of the system is the 4 week battery 
life of the VCO Gelyte batteries. If this modification is 
made, the theoretical system operating life would be limited 

by the six month life of the battery powering the preamplifier 
at the coil. The 14 dB amplifier and the VCO could also be 
integrated into one chassis, thus eliminating one power 
connection and enhancing the reliability of the system. 

A third improvement could be to increase the area and/or 
the number of turns in the coil sensors. This would increase 
the sensitivity of the system, but again, would preclude 
interchange with the underwater system. 

A final area of system improvement could be made in the 
processing of data. Each graph in this research represents 
53-4 man-hours with the majority of time consumed by manually 
reading the data from the spectrum analyzer and then entering 
the data, point by point, into the computer for processing 
and plotting. Since the spectrum analyzer provides the 
data in a digital format, an interface between the analyzer 
and the computer should be possible. Perhaps a Hewlett- 
Packard spectrum analyzer could be obtained with the proper 
frequency range. In this case the interface problem would 


be greatly reduced, since the adapters are already in place. 
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APPENDIX A 


Peewee oe QULPMENT AND TESTS 


A. EQUIPMENT CONFIGURATION 

Most of the equipment utilized in this project was either 
designed and/or manufactured at the Naval Postgraduate 
School, while the test and telemetry equipment is mostly 
of commercial origin. A great deal of the effort expended 
in this thesis was in testing and ensuring that the experi- 
mental equipment operated as conceived and predicted. For 
example, it 1S important that all sensor equipment be dc 
coupled to avoid attenuation of low frequency signals. A 
discussion of the equipment and its configuration follows: 

1. Data Acquisition Equipment 

SENSOR 
The general hookup of the data acquisition was shown 

in Figure 5. The sensor is a self-supporting, continously 
wound, non-center-tapped coil antenna manufactured by Elma 
Engineering, Palo Alto, California, from approximately 5460 
turns of 18 gauge copper magnet wire. — The sensitivity of 
the coil, as calculated and confirmed by measurement, is 
given in Appendix B. Its weight is approximately 50 kg. Its 
overall dimensions (.15m by .37m) are determined by the largest 
glass sphere that 1S commercially available. These spheres are 


used to enclose the coil during underwater experiments. The 
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womlereststance 1S £24 Ohms and its self-inductance is 


BEpouoximately 9.51 henries. 


PREAMPLIFIER 

The differential preamplifier Model 13-10A is an 
improved version of Model 13-10 which was employed by McDevitt 
and Homan (1980). This preamplifier was selected because 
of its low power requirements, compactness and excellent 
low frequency noise characteristics. The final stage of the 
preamplifier also contains an active filter which provides a 
Sharp cutoff for frequencies above 20 Hz. The overall pre- 
Miperraer gain for inputs of less than 2.5 millivolts is 
60 dB. The preamplifier is buried with the sensor coil at 
the test site. It is connected to the rest of the system by 
90 meters of armored cable, for rodent protection. Both ends 
of the cable are protected by diodes to prevent large voltage 
surges due to lightning from entering the system. An 
additional 14 dB amplifier was placed at the end of the cable 
to provide an average voltage range of 0-10 millivolts into 
the VCO. A graph of the overall amplifier gain characteristics 
is included as Figure 29. A schematic diagram of the pre- 


Piper Ter “Circudt 1s “shown in Appendix on 


VOLTAGE CONTROLLED OSCILLATOR 
A schematic diagram of the voltage controlled 
Ssciiiaror {VCO) is included in Appendix C. The function of 


the VCO is to convert low frequency analog voltages received 
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from the sensor through the preamplifier into a time varying 
frequency, centered at 1500 Hz. The output level of the VCO 
iS approximately 0.4 volt (zero to peak) square wave of 

mdmying frequency which is then transmitted by the telemetry 


System. The conversion factor for the VCO is 50 Hz/mV. 


Slit RANoMiSstrTING EQUIPMENT 

An MX Portable 302 series (Motorola Model H23AAU1120N) 
mieeeiannel radio had been selected for utilization at the 
monitoring site, (Beliveau, 1980). It has an output power of 
ewatt, RF bandwidth of 16 KHz, and a channel spacing of 25 
mews lt 1S designed to Operate at 7.5 volts. This radio 
meduired MOdirilcation to operate in the continuous mode for 
telemetry data transmission. A Motorola vehicular antenna 
adapter cable (Part No. NKN-6223) was added to allow the use 
of an external antenna. <A Motorola external microphone 
speaker (Part No. NMN-6070) also was installed. This allowed 
PiewuadioO to be Operated for voice transmission or switched 
to continuous transmission of the VCO output signal. 

A six element Yagi antenna, DB-292 (manufactured by 
Decibel Products, Inc.) Was selected for use with the radio. 
The antenna consists of a reflector, driven element (folded 
dipole), and four directors mounted on a common support boom. 
Gain and VSWR are essentially constant (maximum gain of 9.5 
dB at mid frequency) across a 10 MHz band (providing optimum 
performance if a multi-frequency system is later installed). 


The DB-292 has a beamwidth of 50 degrees. 
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meeieoeseceiving Equipment 


Another Yagi antenna was mounted on the roof of 
Spanagel Hall at the Naval Postgraduate School. The antenna, 
mom, Veetbel Products, Inc., 1s a DB-292-2 (2 stack array) 
eeoviding a 12.1 dB gain. 

A 40 watt output power, two channel, base station 
(Motorola Model L4355B1136) is the transmitter/receiver at 
Spanagel Hall. It has a receiver sensitivity (with preamplifier) 
Sev. 2o Microvolts and an output impedance of 50 ohms. This 
yields a receiver sensitivity of -199 dB. Comparing this 
receiver sensitivity to the minimum monthly median signal 
level expected (as calculated by Beliveau, 1980), yielded a 


Signal excess of about 40 dB. 


PHASE LOCK LOOP 

Radio waves are affected by many factors during their 
travel from transmitter to receiver. Factors such as 
frequency, distance, antenna heights, curvature of the earth, 
atmospheric conditions, and the presence of intervening 
obstacles such as hills, buildings, and trees all enter into 
peewecaleculation of the total loss of the propagation path. 
Each radio system has a maximum allowable transmission loss 
which, if exceeded, results in poor or total loss of reception. 

The Chew's Ridge to Naval Postgraduate School 
transmission is subject to random bursts of static noise 


Superimposed on the audio data tone. These static bursts 
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wees Receiving Equipment 


Another Yagi antenna was mounted on the roof of 
Spanagel Hall at the Naval Postgraduate School. The antenna, 
mise byebecibe! Products, Inc., is a DB-292-2 (2 stack array), 
meoviding a iZ-1 dB gain. 

A 40 watt output power, two channel, base station 
(Motorola Model L4355B1136) is the transmitter/receiver at 
Spanagel Hall. It has a receiver sensitivity (with pre- 
amplifier) of 0.25 microvolts and an input impedance of 50 
Suiiemee 11s yields a receiver sensitivity of -119 dBy. 
Comparing this receiver sensitivity to the minimum monthly 
median signal level expected (as calculated by Beliveau, 1980), 


yielded a signal excess of about 40 dB. 


Pass LOCK LOOP 

Radio waves are affected by many factors during their 
travel from transmitter to receiver. paecore such as 
frequency, distance, antenna heights, curvature of the 
earth, atmospheric conditions, and the presence of inter- 
vening obstacles such as hills, buildings, and trees all 
Tima Geimomuncomeal culation of the total loss of the propaga- 
tion path. Each radio system has a maximum allowable trans- 
mission loss which, if exceeded, results in poor or total 
loss of reception. 

The Chew's Ridge to Naval Postgraduate School transmission 
1s subject to random bursts of static noise superimposed 
On the audio data tone. These static bursts create 
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numerous high frequency components and degrade the system's 
ability to reduce the data to meaningful information. 
Pepndse-lock loop Was designed to track the 1500 + 500 Hz 
audio signal and to filter out sudden random static frequencies. 
The loop parameters are: 

1. Equivalent loop noise band width 500 Hz 


meee damping ratio of 1/ 2. 


TAPE RECORDER 
The tape recorder employed for this project is a 

Hewlett-Packard model 3960 Series Instrumentation Tape 
Recorder. The tape speed was 3 3/4 ips. Additional charac- 
Metistics, £or this tape speed, are: 

imeerreducheyeresponse: S0edze- 15 KHz 

2. Signal-to-noise ratio: 38 dB 

5S. Flutter % peak to peak: 0.4 
Fiutter compensation is incorporated in the recorder to 
compensate for any flutter component introduced by tape 
eeeeaeVariatrons. When the rluttem compensation switch 
located behind the access door on the upper left frame of 
the chassis is set to "On", channel 2 cannot be used for 
fica tTeeenamme. During reproduction, the output of channel 
2 is inverted (phase changed by 180°) and added to the 


other FM channel, thereby cancelling the flutter in these 


channels. 
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Other systems in use at the Naval Postgraduate School 
employed portable cassette type tape recorders that are 
operated while sealed in glass spheres at depths of up to 
300 feet below the ocean surface. To reduce WOW and flutter 
in these recorders, a 2 KHz reference oscillator signal is 
recorded simultaneously with the VCO data signal, and the 
difference frequency is then the input to the frequency-to- 
voltage converter. The reference and data signals are mixed 
on playback. This method was tried but showed no improvement 
relative to results obtained with the HP 3960 with flutter 


compensation. 


FREQUENCY-TO-VOLTAGE CONVERTER (ANADEX PI-375) 
The purpose of the frequency-to-voltage converter 

(FVC) is to convert the frequency difference from the standard 
mgdeHz VCO output into an analog voltage. The varying 
voltage is then a representation of the geomagnetic fluctua- 
tions detected by the sensor. The FVC is powered by 24 Vdc 
and has the following additional characteristics: 

1. Input frequency range: 500-5000 Hz 

meeoOueput de voltage: 0-8 volts full scale 

SeemversioOneractor: 5 mV/Hz 
The conversion factor depends upon the amount of dc offset 
introduced by the FVC. The greater the offset, the greater 
the Mv/Hz sensitivity. The sensitivity possible is limited 


by the need to keep the response linear over the frequency 


68 








range expected. The output of the FVC is then analyzed with 
meopeectrumeanalyzer,. A graph Of the VCO-FVC gain characteristics 


ms Figure 30. 


SPECTRUM ANALYZER (S/A) 

A Mini-Ubiquitous Model 400 was employed for spectral 
analysis. It allows for de coupling, has a digital display 
of the power spectrum, and has a logarithmic, as well as a 
linear frequency scale. 

STwecweiesalSeriminatoOr Output is a positive dc 
voltage, it is necessary to provide for a de offset in order 
to increase the input sensitivity of the spectrum analyzer. 

A Textronic Model AM 502 differential amplifier was used to 


remove the dc offset. 


BO Wieo UO lL ES 
1. Batteries 

Gould Gelyte batteries have been utilized to meet 
the dc power requirements of the preamplifiers and the voltage 
controlled oscillator. They are sealed, rechargeable, and 
virtually maintenance free batteries, containing no free 
electrolyte. The preamplifier is powered by two 12 Vdc 18 
ipenc hour batcertes providing + Vdc. The 14 dB amplifiers 
and the VCO are powered by two 6 Vde 18 ampere-hour batteries 
Fueovtaame + Omvde. These batteries require replacement and 


recharging approximately every three months. 
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The MX 320 radio is powered by 4 low discharge lead- 
acid batteries (Willard Type DH-5-1) connected in series to 
provide 8 volts. Each battery provides 600 amp-hours at 77° 
at 2 volts. Since the geomagnetic data is transmitted 
continuously (100 % duty cycle), the battery system is designed 
to allow for about 28 days operation without recharging. 

This is based upon a power requirement of 22.1 amp-hours per 

day (19.2 amp-hours with a 15% safety factor). The batteries 

would then be completely discharged and require replacement. 
ee OO LameAnray 

The purpose of the solar array is to provide an 8 volt 
source to recharge the continuously discharging batteries. 
Since the array output will vary depending upon the sunlight 
available, an 8 volt regulator is utilized to prevent battery 
damage due to charging voltages higher than 8 volts. This 
8.0 volt rating is a nominal voltage rating and the actual 
voltage output will vary from 9 to 11 Vdc depending on the 
thermistor temperature. The thermistor is a remote assembly 
that senses the ambient temperature in the vicinity of the 
storage batteries and determines the desired output voltage 
Oeeene regulation. 

In addition, a separate low voltage Schottky diode is 
incorporated in the regulator to prevent battery discharge 
Enneuene the Solar array during non-sunshine periods. The 
solar array/battery arrangement is designed for continuous 


operation on the basis of 8 sunless days. The absence of 
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sunlight beyond the 8 day limit will lead to system shut-down 
and battery replacement or recharging. 

The solar panels (MSP43A40, Motorola Solar Systems) 
consist of a series of photo-voltaic cells each providing 
a specific voltage per cell depending upon the amount of 
sunlight available for the conversion of light energy to 
electrical energy. The array consists of two modules in 
parallel. The array was faced toward true south at a slope 
of 47 degrees from the horizontal in order to capture the 
maximum amount of sunlight as the earth revoles and the sun 
"moves'' from the eastern to the western horizon. The amount 
of total energy can be greater if the angle is changed to 61° 
in the winter and to approximately 45° in the summer. In 
general, the months of November, December, and January will 
have the lowest system output to load ratio since the 
temperature will be seasonally low creating a higher rate of 
discharge (poor battery efficiency). The daylight is also 
shorter during these months, resulting in less sunlight 
captured by the solar array. Table 1 shows the predicted 


otemeperrOrmance Lor the various months of the year. 
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APPENDIX B 


SYSTEM CALIBRATION 


Pee OENSOR SENSITIVITY 
1. Theoretical Sensor Sensitivity 
Faraday's Law of Induction shows that the electromotive 
force (emf) induced in a single-turn circuit is equal to minus 
the rate of change of flux. 


-99 


ot. 


Where the flux passing through a surface bounded by the 
conductor is defined as: 
> = fB- Ga 


For N turns of wire connected in series, the total 
emf induced is the sum of the emf induced in each separate loop: 


a® 


emf = “NG 


Assuming that magnetic fields vary sinusoidally with 
time and may be written as: 

B= Bie Jot 
Where’ B| iS a Magnetic vector field that is only a function 
of spatial coordinates and not of time. Therefore we have: 


Hie welsmelcudnea enetosed by the coil, and B, is the axial 


component of B,. 


“J 
GA 





The magnitude of emf is 
Steeler (at)sA 


Since the output of the sensor is connected to a high 
input impedance amplifier, then negligible current can flow 
through the sensor which might otherwise produce a back emf 
and reduce the magnitude of the field to be measured. The 
sensor's open circuit voltage or sensitivity can be calculated 
for the coil parameters given in Appendix A for a InT field 
at a frequency of 1Hz; emf = 2.8316 yV, therefore 20 log(emf) 
Soe. 95 dB(re: 1 volt). 

The area used for the calculation was the average of 
the areas found by using the inner radius and outer radius 
ere the coil. 

2. Experimental Determination of Sensor Sensitivity 

As was found by McDevitt and Homan, (1980), the need 
for a uniform magnetic test field required the use of a .6l 
meter radius Helmholtz coil. The Helmholtz coil produced a 
known field at a specified distance and the electromotive 
force produced by the sensor was then measured. The coil 1s 
designed to produce an easily accessible, almost homogeneous 
magnetic field in the central region between the coils. 

Ewumecsmons stor Ene Magnetic induction vector for 


the axial and radial components are: 
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as given in (Jefemenko, 1966). 


Therefore, in the central region 


ae SNI 
Z SV5 RB 


to the terms of the order (L/R) where L is the linear 
dimension of the region. Specifying the desired field 
strength and distance determines the current through the 
coil at different frequencies. To produce the required 
current, a variable resistor was placed in series with the 


coil and the signal voltage was set according to: 


Nigewe nk, = Coll impedance 


Ro 


variable resistor 

The current was determined to be within + 1.0% of one 
mA and the field produced along the central axis of symmetry 
was determined to within + 10% of a nT. The main reason for 
the slight difference in sensor sensitivity obtained from 
that of McDevitt and Homan was the use of N = 5460 turns 
versus N = 6000 turns, where the first was the actual count 


obtained from the manufacturer. 


is 





Pee REAMIPEIPIERS, VCO, AND FVC 

Once the sensor sensitivity had been verified using the 
Helmholtz coils, the other gains and conversions of the 
electronic component were determined. This was done by 
mipuecting @ known signal at varying frequencies to each 
Somponent and measuring the output. 

ihe Major problem encountered was producing the stable 
calibrated known signal. The source used for the final 
calibration was a Hewlett-Packard 3320B Frequency Synthesizer. 
The main reason for its selection was the output level control 
memes Of +26.990dB to -69.99dB. This allows the input to the 
preamplifier to be in the 1-2 milli volt range. 

The outputs are then measured by oscilloscopes and the 


gains calculated. 


fe) SPECTRUM ANALYZER 


The spectrum analyzer, when in the dBV display mode read: 
N(dBV) = 20 log A volts 
JL sive 
Since 0 dB referenced to InT not 1 volt is desired: 
Revolt s Nevoles/n1 E A volts Mevolts 
wou tC K WOlts/nt Re VOuts/ne ievyouet ni 
Where B (the field measured by the svstem) is: 


oe aioe 
A WGwesS on! 
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so referenced to InT we have: 


B= 


N(dB) = 20 log tnt 


+ 20 log X 


where X is the system transfer function given in units of 
wolts/nt. 

Since spectrum analyzers use various averaging windows 
and have various numbers of bins, the following corrections 


must be applied: 


B? = 6(BW). 


number of bins 


and (BW) = ae) (window corrections) 


Therefore 


ie ie 
N(dB) = 20 log [67?(BW) 7] + 20 log X 


or 
10 10g $= N(dB) - 20 log X - 10 log (BW) 


where 10 log $6 is the corrected reading or the power spectral 
density. 6 has the units of nT /Hz. 

For the Mini-Ubiquitous Model 440, the window correction 
1s 1.5 dB and it uses 400 averaging bins. The BW correction 


on the 0-20 Hz scale was: 
: 20 2 
BW = (z50—) (lS), =" 0075 


tt loco, = Ll.Zs5 
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APPENDIX C 


BOULPMENT SCHEMATICS 


Figure 33 shows the preamplifiers utilized in this 
project. It is an improved low noise ELF amplifier model 
imo UA designed to provide the quietest operation, i.e., 
Mmmeeatest Sigmal-to-moise ratio. Figure 34 is the VCO unit 


meed in this system. 
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APPENDEX D 


EQUIPMENT USAGE AND MAINTENANCE 


The equipment utilized in this sytem requires a certain 
amount of periodic maintenance and handling. Below is a 
list of hints and suggested activities which might help to 
prevent equipment damage and prolong the life of the system. 
PeeeeSeEMOTE SITE EQUIPMENT 

The first action upon arrival at the remote site should 
be to open the storage cabinet and inspect the radio, VCO, 
storage batteries and solar system. (Note: Opening the 
storage cabinet requires a 3/8 inch socket to remove counter 
sunk lag bolts holding the access panel.) 

1. Willard Charge Retaining Batteries Type DH-5-1 require 
the addition of water once a year under normal operations. 
If the solar charging panel voltage is providing more than 
9.2 volts, gassing may occur and more frequent water addition 
may be required. The electrolyte level should be 1 3/4 
inches (4.4 cm.) above the battery plates. The state of 
charge may be checked by using a hydrometer to measure the 
specific gravity of the electrolyte. The hydrometer used 
Should have a float calibrated to read lower than the 
standard automotive type hydrometer. A fully charged cell 
will have a specific gravity reading of 1.300, while a 
discharged cell will read 1.075. The cut-off voltage is 
1.75 volts per cell or a system voltage of 7 volts. 
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The storage batteries supplied with the system contain 
Peamnomeowacic electrolyte. Ihe following general precautions 
for the safe handling of the storage batteries should be 
observed: 

a. In the operation of an acid battery, hydrogen 
gas is liberated which may be explosive if ignited. The 
battery area should be checked for signs of hydrogen gas 
build-up (particularly at the ceiling) prior to performing 
any electrical checks that could be spark-producing. 

b. Plastic or rubber gloves, an apron and safety 
Poeeles should be worn 1f the batteries require handling. 
Fresh water should be available so that the electrolyte 
can be rinsed off immediately if it comes in contact with 
the skin. A mixture of 1 pound bicarbonate of soda in one 
gallon of water can be used to neutralize acid spills. Apply 
the solution until bubbling stops, then thoroughly rinse 
the area with water. 

c. When diluting an electrolyte, always add the 
ueeerglytemto water Slowly, with constant stirring. Explosive 
reactions can occur if water is added to the electrolyte. 

2. Solar Systems 

It is recommended that period maintenance on the solar 
system be performed annually during the springtime. The 
advantage of springtime maintenance lies in its benefit to 
the storage batteries. If water additions are required, 


the battery electrolyte level will be topped off just before 





Pie Tiiomenemseason Of Sreatest water usage - summer. As 
Water is consumed, the proportional acidity of the electro- 
mare wlevel Will inerease and will a achieve its highest level 
during the winter months, providing additional freeze 
mpeotection. 

The minimum yearly required solar array maintenance is 
listed below: 

a. Inspect the solar modules for dirt, dust, bird 
meemmings, Imsects Or other toreign material on the glass. 
Clean with warm water and a soft cloth as necessary. 

b. Inspect the module mounting and anchor bolts 
for tightness and tighten as necessary. 

c. Inspect the silicone weatherization and leads 
of each module for adequate weatherization. Tighten the 
terminals and repair the silicone weatherization as necessary. 

d. Inspect the wiring isulation for cracking, 
peeling or other defects. Replace as necessary. 

Se mcmectetne Installation Site f¢r any new growth 
of weeds, shrubs, trees, etc., which have obscured the 
exposure of the solar array to the sun and which may have 
caused localized shading of the array. Inspect the voltage 
regulator tor bird or rodent nests also. All such obstructions 
should be removed or cleared as necessary. 

f. Note the ambient temperature in the vicinity of 
the thermistor housing. Using a voltmeter capable of 2% 


accuracy at the measured vcltage, check the voltage between 
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the positive battery terminal and the negative battery 
feomminadlweethe Heading Should be between 9.1 and 10.3 volts dc. 

Par eeeincmVvoltace 1S lower than 9.1, the batteries 
should be disconnected and tle voltage measured again. It 
should now be within the prescribed limits. (Explanation: 

If the battery is fully charged and the solar array is 
generating its full output, the regulator will be operating 
to clamp the battery voltage at tke level required for the 
Deeercular battery temperature. If the battery is accepting 
a substantial charge current however, the output voltage will 
be lower and it will be necessary to disconnect the battery 
to force the regulator into operation. 

i piteemes voltage 1s heh, or if it is low even*with 
the battery disconnected, the thermistor ccntrol needs 
fteaaaguicinent. lf this 15 the case, the regulator should be 
Temerned to the factory for adjustment. 

1. Next, measure the short Circuit current of the 
array using a current meter with a sufficiently low internal 
resistance (i.e., its terminal voltage). Disconnect the 
array from the battery and connect the ammeter across the two 
levels. When (scone ching sue mMeber aiter the reading is 
taken, the connection should tbe broken rapidly since there 
will be a tendency to draw an arc and burn the leads. The 
short circuit current reading should Le greater than 4.0 


amperes. 


DD 
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oe eeame enna 

Me@imeecmance Orethe £16¢ld antenna consists of checking 
een wedtherization Of the cable connections between the radio 
and the antenna, insuring that the anchoring system is tight 
and verifying that the antenna is properly aligned with the 
base station. 

Toe leetronic Equipment 

Mitmitiemrculdirecencemmewith) the electronic equipment 
is the proper application of the power supply voltages. 
Improperly applied voltages may result in integrated circuit 
emey7Or printed circuit board damage. The most likely time 
Moameiis to Occur 15 during the replacement of the Gelyte 
batteries that power the VCO and Preamplifier. The suggested 
application of dc power is to first, connect the non-common 
ReMelCi¢> ew OCI Genhtain woRObserve the correct polarities; 
then, connect the common lead, being careful not to short 
the common lead against the equipment chassis. CIS LAE 
power should not be applied unless there is a load on the 
‘input. During normal maintenance the following sequence is 
recommended: 

a. Open the equipment storage cabinet and allow 
any build up of battery gas to displace while checks and 
Service Of the Solar array and antenna are performed. 

Bb. Disconnect the VCO from the radio and use head- 


phones to monitor the sensor signal coming from the VCO. If 
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the signal is the characteristic wavering 1500 Hz centered 
tone, the system 1S operating normally and only battery 
replacement and normal housekeeping are required. If, on 
the other hand, the signal is steady, centered well away 
from 1500 Hz or there is no signal, each component from 
the VCO to the coil must be checked for proper operation. 
A Simpson Model 260 volt-ohmmeter or its equivalent is 
helpful in isolating the source of failure. 

c. For normal service, the next step is to dis- 
Semigece tne + 6 volt battery pack from the VCO, then the 
connector from the armored cable. This 1s necessary since 
any tools, watches, etc. that are carried near the sensor 
coil will cause a large signal that could damage the VCO 
if it was still operating. 

d. The sensor site should be approached by following 
the armored cable, since this route avoids the poison oak 
Ciaeeis Ehroughout the area. The sensor container should 
be checked to insure that the top is high enough above 
ground level to prevent water from storms or melting snow 
from entering. After the cover is removed, a check for 
wildlife is recommended prior to attempting work. Once 
the sensor container is cleared and cleaned, the coil 
Smieitderon smoula be veritved, Next, the power cord is 
disconnected from the preamplifier, and the preamplifier 
batteries, which are located in another container buried 


approximately 15 feet south of the sensor, may then be 
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Powe tGec mater utiiizing the voltmeter to ascertain the 
Pome mmOlertey, the preamplirier is reconnected. Heating 
@uctetape er One Of a Similar type should be used to seal 
the sensor container and plug any access holes. 

e. The VCO may now be reconnnected and powered. 
Head phones are then used to verify proper operation. If 
ememvyCO 1s Operating satisfactorily, the VCO may be connected 
fo the radio. 

A complete set of back-up electronic components, 
a rechargeable soldering iron and spare connectors will 
repair most failures. 
5. Gelyte Batteries 

The batteries used to power the VCO and preamplifier 
are capable of numerous rechargings. Each battery is labeled 
and a record of chargings is maintained. The recommended 
charge procedure stipulates that both the .6 and 12 volt 
18 ampere-hour-rated batteries should be charged at 1 ampere 
poueapproximately 20 hours for a full charge. Do not over- 
charge in excess of 24 hours more than the ampere-hour 
ipa te . 

6. Radio Transmitter 

Other than protection from water and dust, the only 
maintenance required on the MX320 radio is to insure the 
proper polarity of the power supply if it is disconnected 
and then reconnected and to attempt to maintain a constant 


voltage input below 8 volts dc. 
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Eee ok, EQUIPMENT 
1. Radio and Antenna 
The Motorola base radio has no operator maintenance 
requirement. The only maintenance check for the receiving 
antenna is to be certain that there is adequate water pro- 
tection at the connection between the antenna lead and the 
coaxial cable running to the base radio. 


Z. Hewlett-Packard 3900 Series Instrumentation Tape 
heconden 


The most important preventive maintenance considerations 
are: preventing oxide build-up in the transport tape path; 
keeping the recorder clean; and demagnetizing head and tape- 
path components. The main procedures are listed. 

a. lape Path Cleaning 

Note: Do not use toluene or xylene cleaning 
agents on the pinchroller and do not allow cleaning agents 
to come in contact with the magnetic tape. Use cleaning 
agents sparingly and keep them out of bearings, especially 
meeenesDase Of the capstan. 

All transport parts that come in contact with 
the tape require cleaning. Oxide build-up can cause degrada- 
tion of high frequency response, distorted record and play- 
back and tape and head damage. Worn tape heads cause changes 
in operational characteristics, too. Record sensitivity 
increases and there is a decrease in high-frequency response 


as the heads wear. The tape passing over the head presents 
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eetcmrvely sroueh suriace €6 the head. Since even high 
Guality tapes will eventually cause wear on the heads, the 
pie sadetense 1S £0 keep the tape-contacting components clean 
and as free of foreign matter as possible. 

Clean tape heads with cotton-tipped applicators 
dampened in a head cleaner like DuPont Freon TF. Rotate 
Mpemapplicator against the head surface and note any dis- 
coloration of the cotton. Use as many applicators as required 
to avoid wiping the head surfaces with a dirty applicator. ' 
When the applicator comes away clean, the head surface is 
clean. Applicators should also be used to clean the tape 
mudes, rollers, capstan and pinchroller. Take extra care 
to clean the undercut edges of the guides and roller grooves. 
Accumulations in the corners of the guides and roller grooves 
Will damage the tape. 

b. Demagnetizing Heads and Tape Path Components 

The heads and tape path components should be 
Gemapmetized after every 72 hours of operation. To 
demagnetize: 

1. Disconnect the power from the recorder. 
Z. Remove any tape on the recorder. 
5. Remove the head cover and interhead shield. 
4. Hold a bulk tape eraser as close as possible to the 


capstan assembly components. 


a0 





Petes cuanagenold the “On” button and slowly rotate the 
bulk tape eraser in an arc over the capstan assembly components. 
Slowly draw the bulk tape eraser outward until at least two 
feet away, then release the "On" button. 


eaeecoldce tie wimtennedad Shicid and head cover. 


wee 





200 
218 
228 
238 
24D 
258 
e60 
270 
280 
298 
308 
318 
328 
338 
342 
3958 


APPENDIX E: COMPUTER PROGRAMS 


A DATA PLOTTING PROG 

AS: TITLE 

Al$: SUBTITLE 

B$: X AXIS 

C$: Y AXIS 

P =: THE PLOTTING PARAMS 
B—-X AXIS VALUES 
1-AMB BINS 
2aeineErPS7CRASS 

| 3-AVE AMB 

| 4-7 CP ASSES 

| NS¥ #80F PIS - 1 

| N8: # OF PLOTS ON GRAPH 

[| NS: # OF PLOTTING PARMETERS 

| Q$: AN ANSWER 

| QO$: LOG PLOT ANSWER Y/N 

! L$: LOG SCALING FLAG 

| T$e TITLE FLAG 

! F1l$: FLAG SIGNALS FIRST RUN 

DmiMisG: L1, Le, H1, He, I, J, xX, Pl 

Pecra, Ke 

| AS$: AUTO MODE FLAG 

| 

| 

| 

| 


Rea@e KiCe>) |! LOG VALUES 

SHORT D¢2Q8,5) ! THE DATA 

DIM A$(l1@@], A1i$l19@gl1, B$ (1941, C$C190], MS8$C25] 
16 RESET" 

CLEAR 

FiS="RESET" 

DISP “NOTE: ALWAYS ANSWER PROMPTS WITH" 
Cie or OR N FOR ‘YES’ OR “NOQ’®* 


OZ 





360 
378 
38d 
390 
400 
41d 
420 
43G 
449 
430 
460 
476 
480 
49d 
30d 
312 
520 
938 
944 
358 
360 
378 
382 
390 
682 
614 
620 
634 
640 
65d 
662 
670 
680 
690 
700 


Peorren 1S 705 

GOSUB 55d |! READ DATA 
GOSUB 854 ! SELECT&REV DATA 
IF P1l=@ THEN 441@ ! END PROG 
IF Pi=6 THEN 2480 ! REDRAW 
IF ASS="SET" THEN 464 


GOSUB 1380 ! INIT PLOT PARM 
GOSUB 1770 ! DRAW AXES 
GOSUB 2824 ! PLOT DATA 


FouOeseon! MORE FEBS? 
GOSUB 499@ ! AUTO MODE 
GOSUB 1778 ! DRAW AXES 
GOSUB 2828 ! PLOT DATA 
GOTO 382 ! RETURN FOR MORE 


! 
| 
! 
| 
| “READ IN DATA 

N8=1 

NS=e 

LsUasSr §6FLOIS/OCRAPA” 

! INPUT N8& 

eto. FARMS PLOTTED" 

! INPUT NS 

ie Naso. oS THEN 6S2 

DISP "TOO MANY, TRY AGAIN" 
GOTO 68d 

DISP “DATA FILENAME" 

INPUT FS 

ASSIGN# 1 TO FS 

N7=N8+1 


READ® 1, 1 
N7*N7-1 ! SUB FOR MULTI READ 


ao 


a 





—— pee. 
—— er 


712 
723 
732 
74GB 
7354 
762 
778 
782 
792 
822 
81a 
822 
832 
842 
852 
864 
872 
88a 
89a 
992 
914 
922 
932 
944 
952 
962 
976 
982 
994 
1882 
1812 
1022 
1232 
1842 
1852 


IF N7=@ THEN 682 
READ# 1; NS 
NS=NS—1 

FOR I=@ TO NQ 
FOR J=@ TO NS 
READS Tes OC, J) 
NEXT J 

NEXT I 

RETURN 


| PSSEECTOAND REVIEW DATA 


CLEAR 

DISP “CHOOSE PARAMETER" 
DIS @-END PLOTTING" 
esr 1-RAW DATA” 
DISP " Can OCESseD 
VoL 3-UNUSED" 

fm BLSe 4—UNUSED" 

[ DESP ~ S-UNUSED" 
DISP “ SS KRepea LOT” 
INPUT P1 

N6=N5+1 


IF Pl"6 THEN 1322 

IF P1<=N6 THEN 1222 
DISP "DATA NOT AVAILABLE" 

DISP “CHOOSE LOWER NUMBERED PARAMETER" 
GOTO 85a 

F1S="SET" 

IF Pi=@ THEN 1342 

PeP 

ASS="RESET" 
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1262 
1372 
1282 
18928 
1182 
1112 
1122 
1132 
1142 
11358 
1162 
1172 
1182 
1198 
1288 
1218 
1229 
1238 
1242 
1250 
1262 
1278 
1282 
1298 
1300 
1312 
1320 
1332 
1342 
1352 
1362 
1378 
1380 
1394 
1492 


L$S="RESET" 
FiS="RESET" 
DiS “REVIEWS 
INPUT Q$ 
IF Q$<>"yY" 
Hi=D ¢@, a) 
H2=D (2, PD 
Li=D0 (2, a) 
L2=D (2, Pd 
FOR I=@ TO NO! FIND MAX/MIN 
Hi=MAX (DCI, 2), H1) 
LI=MINCDCI, @,L1) 
H2=MAX CDCI, P), Hed 
L2eMIN CCI, Pd, Led 
becky CL, @), DCI, Pd 
NEXT I 

DISRSUStING "3/~ 
DISP “XMIN=";Li; " 
DisrPr USING. 3/~ 
DISP "“YMIN="; Le; " 
DisP USING “3s/" 
DISP "CONTINUE" 
INPUT Q$ 

IF Q$<>"Y" THEN 
DISP “AUTOMODE" 
INPUT Q$ 
IF Q$<>"Y" 
AS$S="SET" 


RETURN 
| 


THEN 1382 


XMAX="3 H1 
YMAX="; H2 
872 


THEN 1342 


! INITIALIZE PLOT PARAMS 
DISP “INPUT XMIN, XMAX, YMIN, YMAX" 


25 


1412 
1420 
143d 
144d 
145d 
146d 
1472 
148d 
1490 
1582 
on 
1528 
1538 
154d 
155d 
1562 
1572 
158d 
1592 
1602 
1612 
16220 
1630 
1642 
165d 
1662 
1672 
1682 
1682 
178 
L7G 
1722 
173d 
1742 
1752 


iwi eis Ac, 11, te 


picme ~ AxlS LOG SCALE” 
hse “ENTER Y .GR N”™ 
INPUT Qd&$ 

IF Q0$<>"Y" THEN 1512 
X3=2 

X4=X1 


toe ere nin ¥, Yell Co 

INPUT Y3, Y4 

GOTO 153d 

Picea iNew crv, yaNTVY, XINTC, YINTC” 

INPUT X3, Y3, X4, Y4 

ime SET TREN REMEMBER TITLES 

IF T$="RESET" THEN GOSUB 1680 ! ENTER TITLES 
RETURN 


oe = 


| 

| 

DP ewe fITLE & XfAXIS 
Ouse “Title” 

INPUT A$ 

Olsr “SUBTITLE” 

INPUT A1S 

Scheme AXIS TITLE” 
INPUT BS 

Store AALS TITLE” 
INPUT C$ 

IF N8>1.5 THEN GOSUB 3622 | 
TS="SET" 

RETURN 


BNUeRe peo! TITLES 


== 


| 
| 
| 
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1760 
V7 
1782 
1730 
1822 
1812 
1820 
1830 
1842 
1858 
186d 
1872 
1880 
189d 
1900 
191d 
1922 
1932 
1942 
1952 
1962 
1972 
1980 
1990 
2822 
2812 
20208 
28328 
204 
205d 
260 
2070 
eZ8a 
2BS 
2140 


l 

[| DRAW AXIES 

Beige eck, 110, 25, 85 

Golvze Ss 

! DRAW LINEAR ? 

IF Q0$<>"Y" THEN GOSUB 2589 

! OR DRAW SEMILOG 

IF Q0$="Y" THEN GOSUB 2772 

[ 

er ON EABELS 

Meter Wie) /7 2, (arate (Y2~-Y1)/25, ~-2 
CSoUZe. 4 

LORG §& 

EABEE A> 

LABEL A1$ 
Gettmeecana) 72, 1i— CY2-Y1)/8, -2 
LORG 4 

LABEL BS 

DEG 

LDIR 92 

LORG 6 

LORG 5 

meOte Xi] (X2=X1)78, (Y1i+Y2) 7/2, -2 
LABEL CS 

LDIR @ 

RETURN 


| 
| 
| 
| MULTI PLOTS 

IF N7<>5 THEN 2120 


LINETYPE 4 
GOSUB 2312 


of 





Zn 
2124 
2138 
2142 
2152 
2160 
2172 
218d 
2190 
2200 
2214 
2228 
2238 
2240 
2250 
2268 
22708 
228g 
2298 
2380 
2310 
2320 
2330 
2340 
2358 
23628 
2378 
2382 
2398 
242 
2412 
2428 
2432 
2440 
2450 


GOSUB 722 

IF N7<>4 THEN 2162 
LINETYPE 6 

GOSUB 2312 

GOSUB 722 

IF N7<>3 THEN 22208 
LINETYPE 7 

GOSUB 2312 

GOSUB 722 

IF N7<>2 THEN 2242 
LINETYPE 8 

GOSUB 2312 

GOSUB 722 

IF N7<>1 THEN 2282 
LINETYPE 1 

GOSUB 2312 

GOSUB 682 

IF N8=1 THEN RETURN 
GOSUB 3062 ! THE PLOT TITLES 
RETURN 

{ 

| 

| 

| PLOT DATA 36626 

IF Q9$<>"Y" THEN 2382 
MOVE LGT (CD Cd, @>, Od, Pd 
GOTO 2392 

MOVE OC@Z, @>,O0¢@, Pd 
FOR I=g TO NS 

IF Q0$<>"Y" THEN 2432 
Peon) (Cl, 2), O¢1, P> 
GOTO 2442 

PEG) <i a @ Clare) 
NEXT I 

PENUP 
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iv q4 


s 


ie! 
- 


7 
¢ 
. = 


44 =f. 


> - 


e 


5 


oe 
® 

r 
— 





2462 
2470 
248 
24908 
2380 
2518 
2028 
2934 
2942 
2008 
2068 
2978 
2982 
209 
2620 
26104 
2620 
2632 
2644 
2650 
2662 
2672 
268g 
26598 
2/0 
2/712 
2728 
2738 
2742 
27358 
2762 
2778 
2/78 
2/798 
288g 


PENUP 


RETURN 

| FOR MORE GRAPHS 

GOSUB 1772 ! DRAW AXIES 
GOSUB 2832 

GOTO 380 


| 
| 
| 
| 


| DRAW LINEAR AXIES 
SOALE™K1, Xe, Y1, Ye 

AXES X3, Y3, X4, Y4 

LORG 5S 

GSIZe 3.5 

FOR X™X1 TO Xe STEP X3 
meOT X, Y¥1-CY¥e-Y1)/55 
LABEL USING "K" 3 X 
NEXT X 

LORG 8 

mem yey! 10 Ye STEP Y3 
PLOT X1-«X2-X1) /4@, Y 
LABEL USING "K" ; Y 
NEXT Y 

RETURN 


| 
| 
| 
| 
{| DRAWS A LOG X-AXIS 
eee ae Ee INEAR. YSAXIS 


LORG 6 


IF L$="SET" THEN 2842 


oe 


7 
Sr 


- Ti. i‘ 





2814 
2822 
283d 
2840 
285d 
c86a 
2876 
c88a 
2892 
298 
291 
2920 
2930 
2940 
295d 
2968 
2972 
2988 
2998 
3888 
341d 
3828 
3432 
3942 
30358 
3062 
3072 
3082 
3890 
318d 
3114 
3128 
313d 
3142 
3158 


L$="SET" 

X1*LGT (X1) 

X2=LGT (Xe) 

FOR I=@ TO 8 

K1 C1) =LGT (I +2) 

NEXT I 

S@aLe Xl, Xc, Y1, Ye 
MOVE Xi, Yi 

DRAW X2, Y1 

FOR I=Xi TQ X2 

MOE) Io@ne <Ye-Y1)/355 
EaeeiaeUSING “K" 3; 1@71 
FOR J=@ TO 8 

K2=Ki1 (J) +I 

MOVE Ke, Y1 

DRAW Ke, Y1+(Y2-Y1) /12¢ 
NEXT J 

NEXT I 

VARS Xl, YS 

LORG 8 

roOi~ersm 10 Y2 STEP Y3 
PLOT X1-¢X2-X1) /4@, Y 
LABEL USING "K" ; Y 
NEXT Y 

RETURN 


| 
| 
| 


| BRAW PEOT TITLES 

IF FisS="SET" THEN 3228 
Dis wonrnt) COERES" 
UNG CL, Ce 

C3=Ce 





3160 C8=(X2-XK1)/122 

31706 CS9=(Y2-Y1)/1284 

318@ MOVE C1*C8+X1, Ce*xCS+Yl 
3198 DISP "POSIT OK" 

3286 INPUT Q$ 

3210 IF Q$<>"Y" THEN 31192 
3228 Ce=C3 

3232 MOVE (C1+12) *C8+X1i, C2*xC9+Y1l 
3248 CSIZE 3 

3258 LORG 5 

3268 LABEL MSS 

3276 GOSUB 3572 ! STEP DOWN 
3288 IF N8=5 THEN 332g 
3290 IF N84 THEN 3352 
338@ IF N8=3 THEN 3382 
3310 IF N8=e2 THEN 3412 
3320 LINETYPE 4 

3338 LABEL MS$ 

3348 GOSUB 3492 

33528 LINETYPE 6 

33608 LABEL M4$ 

3378 GOSUB 34392 

3388 LINETYPE 7 

3398 LABEL M3$ 

3408 GOSUB 3492 

341@ LINETYPE 8 

3428 LABEL Mes 

3438 GOSUB 34392 

3440 LINETYPE 1 

34528 LABEL Mi$ 

3468 COSUB 3492 

3476 LINETYPE 1 

3480 RETURN 

3498 | 

35288 | 


= 








3o10 
3522 
3930 
3542 
JIo0 
3568 
3572 
3588 
3592 
3642 
3612 
3622 
3638 
3642 


365d ! 


3668 
3670 
3688 
3698 
37082 
3718 
3728 
3732 
37 4G 
3758 
3760 
3778 
3788 
3798 
3880 
3818 
3822 
3838 
3842 
385d 


[ 
| 


| "ORAW LINETYPE ANB MOVE TO NEXT LABEL START 
MOVE (18+C1) *C8+X1, Ce*xC9+Y1 

DRAW (C1+25) *C8+X1, C2*C9+Y1 

CeeCe-5 

MOVE Ci*C8+X1, C2*C9+Yi 

RETURN 


e 


(i) EMER PLOT TITLES 
IF N8<4.5 THEN 3782 
Oise INPUT TITLE MS" 
MSS="RES = 188 HZ" 

fe "Ne=3. 5 TREN 3730 
Oise Net TITLE M4” 
M4$="RES = 75 HZ" 

IF N8<2.5 THEN 3762 
mor linea 11tLE M3~ 
M3S=""RES = 58 HZ" 

IF N8<1.5 THEN 3792 
DISP "INPUT TITLE Me" 
M2S="RES = 25 HZ" 
Giese INPUT TITLE M1” 
MiS="RES = 1 HZ" 
urs "INPBP TITLE MSs" 
MSS="PRF RESOLUTION" 
RETURN 

| 

| 


at)Z 





3868 
3870 
3882 
3892 
39282 
3918 
3928 
3938 
3942 
3952 
3962 
3972 
3982 
39928 
4822 
4012 
4829 
4832 
4B42 
4852 
4062 
4872 
4082 
4092 
4182 
4112 
4128 
4132 
4142 
4152 
4162 
4178 
4182 
4132 
4208 


| 
iene 15 A TEST STUB 
[week ThE TYPE BOX 
INPUT MS$ 

DISP "ENTER MS$" 
INPUT M5$ 

DISP "ENTER M4$" 
INPUT M43 

DISP "ENTER M3$" 
INPUT M3$ 

DISP "ENTER Mes" 
INPUT M2$ 

BISP “ENTER Mis” 
INPUT M1$ 

SCALE @, 10, 0, 18 
X1=4 

X2=18 

Yi=@ 

Y2=12 

FigS="RESET" 

GOSUB 3112 

GOTO 4262 


| 


{ THE AUTOMODE ROUTINE 


CLEAR 


DISP 
BiSr 
BISe 
i Si 
DSF 
DISP 
OLS 


“WELCOME TO AUTOMODE" 


"AUTOMODE INITIALIZES ALL " 
"PLOTTING PARAMETERS EXCEPT FOR” 
"THE TITLE” 





4210 
4220 
4230 
4240 
4258 
4260 
4270 
4262 
4298 
4380 
4319 
4320 
4330 
4342 
4350 
4362 
4372 
4382 
4390 
4400 
4410 
4420 
4430 
4448 
4430 
4469 
4478 
4488 


4492 
4380 
4512 
4520 


INPUT A$ 

X1i=. 01 

X2=18 

X3=0 

X4=X1 

Y1=-62 

Y2=2 

Y3=12 

Y4=-62 

QOS="Y" | LOG FLAG 
BS="FREQUENCY ¢CHZ) " 

C$="POWER DENSITY ¢dBD " 

IF P<>1 THEN 4352 

Al$=" " 

IF P<>2 THEN 4372 

A1l$="" 

PiSre votes eN) EINe= WHEN PLOTTER READY” 
BEEP 

INPUT Q$ 

RETURN 

ASSIGN# 1 TO * 

Piomeore |” 

BEEP 

Sur 

PRINTER IS 793, 62 

PRINT "&1114@F", "SaSL", "Sa20Qm" 
PRINT °" 

PRINT "“’MIKPLT’-Hustom 14 OCT 81", " 


PRINT 

PLIST 108, 4442 
SOF 

END 





10 | 
28 | 
32 C 
40 R 
32 D 
62 I 
78 C 
80 A 
SO A 
102 
112 
122 
UENC 
132 
148 
15d 
162 
178 
18d 
192 
238 
218 
228 
238 
248 
258 
268 
278 
282 
296 
380 
312 
D¢ 
320 
338 


340 
358 


PRINTER IS 723, 62 

PRINT "&11L" 
LEAR 
EAL D(28d, 3>, T<199) 

ISP "FILENAME? (6 OR LESS CHARACTERS) " 
NPUT FS 
REATE FS, 19 
SSIGN# 1 TO F$ 
SSIGN#¥ 2 TO "TRANS" 
READ#¥ 23 TQ 
ASSIGN# 2 TO * 
DISP "TO GET OUT OF ENTER MODE, ENTER NEGATIVE FREQ 
y" 
DISP 
DISP "ENTER CONSTANT" 

INPUT C 
N9=G 
DISP "ENTER FREQ, PARAM: DATA POINT=", NS 
INPUT E(@),EC1) 

IF —(@)<@ THEN 262 
E (2) =—1#E (1) —-T (NG) -C 


FOR J=@ TO 2 

DC CNY, JD =E CJD 

NEXT J 

NS=N9+1 

GOTO 172 

DISP "DO YOU WANT TO EXAMINE DATA? CY OR ND" 
INPUT A$ 

IF AS="N" THEN 412 

FOR I=@ TO NS 

J=¢Z 

PRINT "DC"; I; "@)="; DCI, Os; " OC"s I: "1)="3 DCI, 133 " 
"sI3; "Q="3D0CI, 2 

NEXT I 


DISP "DO YOU WANT TO CORRECT ANY ENTRIES? CY OR ND" 


INPUT BS 
IF BS="N" THEN 412 


Je 





362 


DISP "TO CHANGE ANY ENTRIES TYPE IN DCI, JD=(NEW ENT 


nme CER), 


370 


Oper Witt PRE PNDEXES I,J OF THE ENTRY YOU WANT TO 


CHANGE. ” 


388 
YOU 
398 
482 
41B 
428 
430 
aad 
450 
460 
478 
482 


DISP “WHEN COMPLETED CHANGING AS MANY ENTRIES THAT 
WANT HIT THE CONTINUE KEY" 

PAUSE 

GOTO 262 


NEXT J 

NEXT I 

ASSIGN# 1 TO * 
END 





18! PRINTER IS 783, 62 

2Oeeeren Sli” 

SOeouMe ll 6199) 

40 CLEAR 

BUPeREATE “1RANS', 19 

62 ASSIGN# 1 TO "TRANS" 

78 DISP "THIS ALLOWS YOU TO ENTER TRANSFER FUNCTION DAT 
ApeeNTO AP BATA FILE. ~ 

soeotse 10 Gel OUT OF ENTER MODE, ENTER 999” 
928 N=2 

108 DISP "ENTER VALUE N=",N 

11@ INPUT T (ND 

1e@ IF T(N)=999 THEN 150 

138 N=N+1 

14@ COTO 142 

158 DISP "DO YOU WANT TQ EXAMINEDATA? (CY OR ND" 
162 INPUT D$ 

170 IF D$="N" THEN 292 

18@ FOR I=@ TO N 

190 PRINT "TC"s3 Is “J="3 TCT) 

200 NEXT I 

218 BiSP “DO YOU WANT TO CORRECT ANY ENTRIES? CY OR ND" 


2249 INPUT C$ 

238 IF C$="N" THEN 290 

248 DISP "TO CHANGE ANY ENTRIES TYPE IN TCID=(CNEW ENTRY 
CER? 

208 DISP " WITH I BEING INDEX OF ENTRY YOU WANT TO CHAN 
Ceow. 

260 DISP "WHEN COMPLETED HIT THE CONTINUE KEY. " 

2/o PAUSE 

e88 GOTO 1528 

290 PRINT# 1: TO 

382 ASSIGN# 1 TO * 

314 END 





BIBLIOGRAPHY 


Barry, J. M., Power Spectra of Geomagnetic Fluctuations 
Between 0.1 and 10 Hz, M.S. Thesis, Naval Postgraduate 


School, Monterey, California 1978. 
Pemevedti fan. lelemetry Design for Solar Powered 


Geomagnetic Monitoring Station, M.S. Thesis, Naval 
Postgraduate School, Monterey, California 1980.. 

Campbell, W. H. and Matsushita, S., Physics of Geomagnetic 
Phenomena, v. 1, 2, Academic Press, 196/. 

Clayton, W. F., Power Spectra of Geomagnetic Fluctuations 
Between 0.4 and 40 Hz, M.S. Thesis, Naval Postgraduate 


Benmooly Monterey, California. 1979. 


Meech onmenwes. GC. and Buxton, J. L., "Superconducting 
Magnetometer Measurements of Geomagnetic Activity in the 


SeemcOnlasizebuequency Range’, Journal of Geophysical 
Peocoamenn ewe oe, Now 27, po. Sl4a=514/, August, T9735. 


Jacobs, J. A., Geomagnetic Micropulsations, Springer, 
New York, 1970. 


Kanamori, K., Takeuchi, H. and Uyeda, S., Debate About 
the Earth, Freeman Cooper and Company, 1970. 


McDevitt, G. R. and Homan, B. B., Low Frequency Geomagnetic 
Pidetiatwons (.04 to 25 Hz) on Land and on the Floor of 
Monterey Bay, M.S. Thesis, Naval Postgraduate School, 


Momeeweye California, 1980. 


Polk, C. and Fitchen, F., "Schumann Resonances of the 
Earth-ionosphere Cavity - Extremely Low Frequency Reception 
Do ieneOnwmn.l. genes. Nat. Bur. Stand. Sect. D., v. 66, 
Dee > See LOZ. 


Salleuugecmmmemen., and Parker, D. C., "The Polarization and 
PoOweteepceenums Of Pe Micropulsations in Bermuda", Journal 


Or Geom icemeammnescarch, v. 68, p. 5545, 1963. 


Wertz, R. and Campbell, W. H., "Integrated Power Spectra 
of Geomagnetic Field Variations with Periods of 0.3-300 s", 


ean OrmmecomlystiGab sReseanch, v. 81, No. 28, p. 3151-5159, 
1966. 


108 





ives tRIBUTLEON LIST 


No. 


Defense Technical Information Center 
Cameron Station 
Alexandria, VA 22314 


Library, Code 0142 
Naval Postgraduate School 
Monterey, CA 93940 


Department Chairman, Code 61 
Department of Physics and Chemistry 
Naval Postgraduate School 

Monterey, CA 93940 


Professor 0. Heinz, Code 61Hz 
Department of Physics and Chemistry 
Naval Postgraduate School 

Monterey, CA 93940 


Professor Paul H. Moose, Code 62Me 
Department of Electrical Engineering 
Naval Postgraduate School 

Monterey, CA 93940 


CPT Michael W. Beard 
300 West Broad Street 
Angola, IN 46703 


Mr. William Audahazy 

Vegieomrp nescarch G& Development Center 
Annapolis Laboratory 

Annapolis, MD 21402 


Chief of Naval Research 
Department of the Navy 
800 North Quincy Street 
AplinegoemeasyA 2221/7 

Code 100C1 

Code 460 

Code 464 

Code 480 


Commanding Officer 

Office of Naval Research Branch Office 
1030 East Green Street 

Pasadena, CA 91106 


Cepaes 


Mm OR 





OM 


PD 


rae 


ee 


14. 


PS 


Dae COT 

Naval Research Laboratory 
Code 2627 

Washington, D.C. 20350 


Office of Research, Development, Test, 
and Evaluation 

Department of the Navy 

Code NOP-987J 

Washanocon, 0.C. 20350 


Paper oree Office Of Sclentific Research 

Department of the Air Force Directorate 
of Physics (MPG) 

Building 410 

Bolling Air Force Base 

taoominmeton, 0.C. 20532 


Army Research Office 

Department of the Army 

Geosciences Division 

Boxed 22.11 

Meccawenmurramole Park. N.C. 2/7709 


PoE Unccmmacerorary Or the Ariny 
fOrpmOpe hat LOlmsekhesearcn 

Moone ce ZO). Pentagon 

Hasnain ton. D.C. 20310 


LTC Larry L. Davis 
Material §& Safety Division 
Ue eNeA, Building 2073 

7500 Backlick Road 
Somimorield, VA 22150 


8 




















Thesis = 

B2895 Beard 

el Power spectra of 
geomagnetic fluctu- 
ations between 0.02 
and 20 Hz. 





